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Abstract 



We use a variety of series expansion methods at both zero and finite temper- 
ature to study an antiferromagnetic Heisenberg spin model proposed recently 
by Miyahara and Ueda for the quasi two-dimensional material SrCu2(B03)2- 
We confirm that this model exhibits a first-order quantum phase transition 
at T = between a gapped dimer phase and a gapless Neel phase when the 
ratio x = J' /J of nearest and next-nearest neighbour interactions is varied, 
and locate the transition at x c = 0.691(6). Using longer series we are able to 
give more accurate estimates of the model parameters by fitting to the high 
temperature susceptibility data. 

PACS numbers: 75.10.Jm., 75.40.Gb 
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I. INTRODUCTION 



The discovery of high-temperature superconductivity has stimulated an enormous 
amount of activity in the study of two-dimensional antiferromagnetism, which may be con- 
nected with the superconductivity phenomenon. The pseudo spin-gap behaviour observed in 
the high T c cuprates has also stimulated intense interest in systems with spin gaps. Several 
new spin gap systems have been found experimentally. Among them some of the compounds 
which have two dimensional character include the coupled spin ladder systems, SrCu203 [|1], 
CaV 2 5 |H, (V0 2 )P 2 7 §, Cu 2 (C 5 H 12 N2)2Cl4 !, and the plaquette RVB system, CaV 4 9 
§■ 

Recently a new two dimensional spin gap system SrCu2(B03)2 has been found by 
Kageyama et al. ||. It has a spin-singlet ground state with a finite spin gap ~ 20 K. 
They also found that the peak of the susceptibility is much suppressed compared with stan- 
dard dimer models, and observed two plateaus in the magnetization at 1/4 and 1/8 of the 
full moment. 

Miyahara and Ueda |7| showed that these observations could be understood on the basis 
of a simple two-dimensional Heisenberg antiferromagnet model with nearest-neighbor and 
next-nearest-neighbor couplings. The copper ions in the SrCu2(B03)2 compound are all 
located at crystallographically equivalent sites, forming a distinctive pattern. Miyahara and 
Ueda J?| show that, remarkably enough, a singlet dimer state forms an exact eigenstate 
of this Hamiltonian at all couplings, and is the ground-state in a region where the nearest- 
neighbour coupling dominates. There is a first order phase transition to a Neel ordered state 
at J' / J = 0.7 ± 0.01. They find that the SrCu2(BOs)2 system lies close to this transition, 
which explains the unusual temperature dependence of the magnetization. The plateaus 
observed in the magnetization curve can also be understood, on the basis that the triplet 
excitations from the ground state are almost localized. Their conclusions were reached on 
the basis of exact diagonalization calculations for lattices of up to 20 sites, and low order 
dimer expansions in (J'/ J), and high temperature expansions in J/T and J'/T. 

In this paper, we reinforce these conclusions by carrying out a more extensive series 
study of the model. These include high-temperature expansions, Ising expansions at zero 
temperature, and dimer expansions both at zero and finite temperature. The model is 
presented in Section II, and the framework for the various series expansions is outlined. The 
results are presented in Section III. It is shown that the temperature dependence of the 
susceptibility can be fitted accurately and in detail by this method. 

II. SERIES EXPANSIONS 

The magnetic properties of SrCu2(BOs)2 may be described by the two-dimensional spin 
S = 1/2 Heisenberg antiferromagnetic model with nearest-neighbor (n.n.) and next-nearest- 
neighbor (n.n.n.) interactions 

H = J y £S i -S j + J' £ Si-Sj . (1) 

n.n. n.n.n. 

The system is illustrated in Fig. 1(a). We denote the ratio of the couplings as x, that is, 
x = J' / J. In the present paper, we study only the case of antiferromagnetic coupling, where 
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both J and J' are positive. In the large J'/ J limit, the model is topologically equivalent to 
the two dimensional nearest-neighbor square lattice Heisenberg model. In the small J' / J 
limit, every pair of spins along nearest-neighbor bonds interact only weakly with each other, 
and the dominant configuration in the ground state is the product state with the spins along 
nearest-neighbor bonds forming a spin singlet. In fact, it has been proved [|8],|7j that this 
state is an eigenstate of the system for any x, and it is the ground state for small enough 
values of x. We can also prove that this system still has this perfect dimer state as an 
eigenstate even if we include the coupling J 2 denoted by dotted lines shown in Fig. 1(b). 

We have studied this system by using various linked-cluster expansion methods including 
dimer expansions at both zero-temperature and finite-temperature, Ising expansions at zero- 
temperature, and high temperature expansions. The linked-cluster expansion method has 
been previously reviewed in several articles [|9| [3~3[] , and will not be repeated here. Here 
we only summarize the expansion methods used, and the results derived from them are 
presented in the next Section. 



A. Dimer Expansions at T = 

At temperature T = 0, we can construct an expansion in x by taking the first term of 
H as the unperturbed Hamiltonian and the second term in if as a perturbation. That is, 
the Hamiltonian of Eq. (0) can be rewritten as 

H = H + xV , (2) 

where 

Ho = Si ■ Sj , 

n.n. 

(3) 

V = £ s 4 • S i , 

n.n.n. 

and where we have set J = 1 for convenience. The unperturbed ground state is then a 
product state of nearest neighbor singlet dimers and the perturbation couples these among 
themselves and with the pair triplet states. As mentioned above, the unperturbed ground 
state is also an eigenstate of the full Hamiltonian, but is not the true ground state for x > x c . 

Dimer expansions can be developed for all ground state properties as well as for the triplet 
excitation spectrum. Here, because of the trivial nature of the ground state in the dimer 
phase, we concentrate on the lowest triplet excitations. We have calculated the dispersion 
relation A(k x , k y ) to order x 15 , extending the calculation of Miyahara and Ueda by 11 terms. 
This calculation involves 11586 linked clusters up to 8 sites. The resulting series coefficients 
are given in Table I. 



B. Ising Expansions at T = 

In the limit that J' ^> J, the model is topologically equivalent to the two dimensional 
square lattice Heisenberg model, so we expect that the system has Need order: an order in 
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which every pair of spins along a horizontal nearest-neighbor bond (denoted as A) has spin 
up, while every pair of spins along a vertical nearest-neighbor bond (denoted as B) has spin 
down. 

To construct a T = expansion about the Ising limit for this system, one has to introduce 
an anisotropy parameter A, and write the Hamiltonian for the Heisenberg-Ising model as 

H = H' Q + XV , (4) 

where 



i£B i€A 



v = 52(Sfs* + sfs]) + * £ (s*s* + sfs]) - * e s- - £ s<\ , (5) 

n.n. n.n.n. igB ieA 

The last term in both H' Q and V is a local staggered field term, which can be included 
to improve convergence. The limits A = and A = 1 correspond to the Ising model and 
the isotropic Heisenberg model, respectively. The operator H' Q is taken as the unperturbed 
Hamiltonian, with the unperturbed ground state being the usual Neel state. The operator 
V is treated as a perturbation. It flips a pair of spins on neighboring sites. 

Ising series have been calculated for the ground state energy per site, E /N and the 
staggered magnetization M, for several ratios of couplings x and (simultaneously) for several 
values of t up to order A 9 . The resulting series for x = 0.7, 0.8, 1, 2 and t = are listed in 
Table II; the series for other value of x are available upon request. 

At the next stage of the analysis, we try to extrapolate the series to the isotropic point 
(A = 1) for those values of the exchange coupling parameters which lie within the Neel- 
ordered phase at A = 1. For this purpose, we first transform the series to a new variable 

5 = i-(i-\y/\ (e) 

to remove the singularity at A = 1 predicted by the spin-wave theory. This was first proposed 
by Huse fl4[] and was also used in our earlier work on the square lattice case fl5fl . We 



then use both integrated first-order inhomogeneous differential approximants [HJ and Pade 
approximants to extrapolate the series to the isotropic point 5 = 1 (A = 1). The results of 
the Ising expansions will be presented in the next section. 

C. High-temperature series expansions 

We now turn to the finite-T thermodynamic properties. We have developed high- 
temperature series expansions (3 = l/(ksT) for the uniform magnetic susceptibility x{T) 
and the specific heat C(T), for the system with arbitrary x, 



. . 1 x - * - Tr5*j Sj c ^ 
1 XU ) = -jj 2^ 1^ Tre-PH 

* 3 

(7) 
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where N is the number of sites, and the internal energy U is defined by 

TT TrHe~P H 

The series for x(T) and C(T) were computed to order j3 7 for arbitrary x. The expansions 
have the following general form: 



OO I 

E[E 

i=0 j=0 



P l , (9) 



where the numerical coefficients. The series coefficients are not presented here, 

but are available upon request. Some of the coefficients can be obtained by using the 
dimer expansion series at finite temperature discussed in the next subsection. The high 
temperature series expansions were first carried out by Miyahara and Ueda |7[ up to order 
f3 2 . Our results agree with these previous results, and extend the series by 5 terms. 



D. Dimer expansions at finite temperature 



The study of finite temperature properties via series expansions is usually done by the 
high temperature expansion method as mentioned above, where we expand in powers of (3 for 
given ratio of exchange couplings. This method often performs poorly at low temperatures 
for many parameter regimes of interest. To overcome this difficulty, one can develop the 
dimer expansion at finite temperature, where we expand the thermodynamic quantities in 
power of x for arbitrary temperatures. This method has been used by one of us previously 
for the Hubbard model |]T2[ and more recently also by Elstner and Singh [TB[] for spin models. 
This method has shown excellent convergence for a wide range of coupling constants at all 



temperatures for the bilayer Heisenberg model |13| and for alternating spin-chains and spin 
ladders 



m 



To get an expansion for thermodynamic quantities, such as the susceptibility and the 
specific heat, in x at arbitrary temperature for the Hamiltonian in Eq.(fJ), one basically 
needs to expand 



e -f3(H +xV) 

in powers of x. This can be obtained by using the following relation 

e -P(»°+xV) = e -^oJ2(-x) n I n , 



(10) 



(11) 



where I n are the integrals given by 

In 

with 



P rti rtn-i _ 

dh / dt 2 ■ ■ ■ / dtnVitJVfo) ■ ■ ■ V(t n ) 
o Jo Jo 



(12) 



V(t) 



D tH 



Ve 



-Wo 



(13) 
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The integrations needed in this expansion are of the type: 

t k e lt , (14) 

and can easily be integrated analytically. 

The expansions for the susceptibility x an d the logarithm of the partition function In Z 
per site take the following general form ]T5| : 



n =o 12n! V12 

(15) 

lnZ(x,/3) = 3/3/4 + lnZ + 5:^M(^) n , 

where Z = 1 + 3exp(— f3), and the coefficient f n (P) for both \ an d hiZ have the following 
general form 

n n+1 

/n(/^) = EE^%- Z , (16) 
k=0 1=0 

where c^j are expansion coefficients. The results up to order n = 6 are given in Table III. 
From this expansion, one can recover the results of the high-temperature expansion and 
dimer expansions if we reexpand in powers of /3, or x at T = 0. 



III. RESULTS 

Having obtained the series for the various expansions above we present in this section 
the results of series analysis. We use integrated first-order inhomogeneous differential ap- 



proximants [16] and Pade approximants to extrapolate the series. 



A. Phase diagram 

The ground state energy per site Eq/N is shown in Figure [| The full points in the large 
J'/ J region are obtained from the Ising expansion, and the horizontal line corresponds to the 
eigenenergy of the perfect dimer state (which is E = 3NJ/8 exactly). These curves cross 
at a transition point x c , which corresponds to a first order ground state phase transition, 
resulting from a level crossing. The numerical estimate of x c = 0.691(6) is a more precise 
estimate of the result x c = 0.70(1) of Miyahara and Ueda discussed above. 

The staggered magnetization M for those values of the exchange coupling parameters 
which lie within the Neel-ordered phase at A = 1 is shown in Fig. 0. We can see that M 
decreases as we turn on J, and appears to vanish at around J/ J' ~ 1.4, near the transition 
point determined from the ground state energies. The errors are too large to determine 
whether or not where is a finite discontinuity at the transition. 
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B. Triplet excitation spectrum 



From the dimer expansions, one can estimate the triplet excitation spectrum for those 
values of the exchange coupling parameters which lie within the dimer phase. The triplet 
excitation spectra for x = 0.6, 0.65 and 0.678 are shown in Fig. f|. We can see that the 
minimum energy gap is at (k x , k y ) = (0, 0) (and the equivalent point (ir, ir)). The bandwidth 
is quite small, indicating that the triplet excitations are highly localized, but increases as x c 
is approached. 

Fig. |5| shows the triplet excitation gap A = A(0, 0) as a function of x. At J 1 = 0, A/ J 
is equal to 1 exactly, corresponding to a single dimer excited to a triplet state. At the first 
order transition point x c , the gap is A = 0.14(5), where the uncertainty is largely associated 
with the uncertainty of x c . 

C. Thermodynamic properties 

Let us now discuss the thermodynamic properties. Fig. |] shows the internal energy per 
site E/NJ vs. temperature /c^T/J obtained from the high temperature expansion and from 
the dimer expansion at finite temperature for x = 0.5. We can see that the direct sum to the 
dimer expansion series at finite temperature converges extremely well down to T = 0, and 
it can recover the exact ground state energy E /NJ = — 3/8 at T = 0. The results obtained 
from the integrated differential approximants to the high-temperature series expansion only 
converge well down to ksT/J ~ 0.3. 

The results for the specific heat at x = 0.4 obtained from the high-temperature expansion 
and from the dimer expansion at finite temperature are shown at Fig. [7]. We choose x = 0.4 
rather than x = 0.5 because the specific heat series converges poorly for larger x. 

Finally, the results for the susceptibility \ at x = 0.5 are shown in Fig. ^. Again, we can 
see that the results from a direct sum of the dimer expansion series at finite temperature 
converges well down to T = 0, (it converges less well near the peak position), while the 
integrated differential approximants to the high-temperature expansion series only converge 
to the position of the peak. 

It is clear that the dimer expansion at finite temperature gives much better results than 
the high-temperature expansion for those values of the exchange coupling parameters which 
lie within the dimer phase. It is interesting to explore this also for ratios of the couplings 
lying within the Neel-ordered phase at T = (x > 0.691). In Fig. |9| and Fig. [l^, we 
present the results for internal energy and susceptibility for x = 1.5, well outside the dimer 
phase at T = 0, and one can see that dimer expansions at finite temperature still give good 
convergence in the high temperature region. Evidently at high temperature the system is 
highly disordered and both methods include contributions from all states. 

D. Compound SrCu 2 (B0 3 ) 2 

Finally we compare the experimental results for the compound SrCu2 (603)2 with the 
theoretical calculations to get an estimate of the exchange constants. Kageyama et al. § 
determined the excitation gap A = 25 K from the NMR relaxation rate. They also measure 
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the temperature dependence of the magnetic susceptibility, and this implies a slightly smaller 
gap: A = (19 ± 1) K. Since there is some uncertainty in the energy gap, to determine the 
exchange constants we begin with the magnetic susceptibility. To convert the theoretical 
x(T) into units emu/Cu mol in experiment, we multiply our \ (calculated with J = 1) 
by NAg 2 fi 2 B/kB, with fis the Bohr magneton, fee the Boltzmann constant, Na Avogadro's 
number, and g the effective ^-factor. There are three fitting parameters: J and x (or J') 
and g. Our goal is try to find a proper parameter set (J,x,g) which gives the minimum 
value of 

P = T,[x exp m)-X thc °(T i )] 2 , (17) 

where x exp ? X theo (^i) are the experimental and theoretical susceptibilities, respectively, and 
the summation is over all experimental points Tj. As the theoretical susceptibility from 
the 6th order dimer expansion at finite temperature converges well above temperatures 
corresponding to a physical temperature To ~ 100 K, we restrict the minimisation function 
to temperatures T > T . We then scan the values of P around the expected parameter region 
for SrCu 2 (B0 3 ) 2 to locate the minimum P min = 1.14 x 10~ 9 (emu/CuMol) 2 at J = 92.0 K, 
x = 0.563 and g = 2.104. This is shown as the full point in Fig. [11] 

Since there is some error in both the experimental and theoretical \i we cannot expect 
that the true parameters for SrCu2 (603)2 are located exactly at the minimum point of P, 
(in fact, the location of the minimum point depends on the chosen value of T ), we need 
to consider the energy gap also. In the 3-dimensional parameter space (J,x,g), there is a 
long tube-like region which has P < 5P m i n . The intersection with the plane g = 2.108 is 
shown as the dashed lines in Fig. 0(a). If we fix the value of J and vary x and g, we obtain 
curves x(J) and g(J) which give near minimum values of P. These curves x(J) and g(J) 
and the corresponding P are given in Fig. II]. We can see that along x(J), P only changes 



slightly by about 10%, and g is about 2.11 almost independent of J. To further determine 
the parameters for SrCu 2 (B0 3 )2, we need to consider the energy gap (which relates to the 
low-temperature behavior of x(T))- From Fig. [5], we can get two curves in (J, x) space 
where the energy gap is 19 K and 25 K, respectively, as shown in Fig. [TTj. We see that they 
cross with the optimal x( J) curve obtained above at 

J = 82.0K, x = 0.678, g = 2.108 

if the energy gap is 19 K, or at 

J = 83.2K, x = 0.664, # = 2.108 

if the energy gap is 25 K. We believe this should be the best fit to both x(T) an d the energy 
gap. 

If we take the energy gap to be 19 K, the comparison of the experimental data with 



theoretical calculations is shown in Fig. [12] and Fig. |13|. We can see the fit is extremely 
good in the high temperature region. We note that the dimer expansion sums tend to show 
sharp peaks, which are not visible in either the experimental data or the high-temperature 
expansions. We suspect that these are an artefact, as the results will become more sensitive 
to the truncation of the dimer expansion at finite order near the transition point in x. The 
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ratio of couplings x obtained here is similar to that obtained by Miyahara and Ueda j7|, but 
they found J = 100 K which is nearly 20% larger than our estimate. They obtained their 
results by considering the paramagnetic Curie- Weiss constant 9. The Curie- Weiss constant 
given by the high-temperature expansions is 9 = ( J + 4J')/4. From experimental data, 
the susceptibility at high temperature can be fitted with a Curie- Weiss constant 9 = 92.5 
K and an effective (/-factor g = 2.14. || The result curve for (J + 4J')/4 = 92.5 K is 
shown as a dotted line in Fig. 11(a). The discrepancy between the two estimates is perhaps 
due to uncertainty in extracting the Curie- Weiss constant from experimental data at finite 
temperature. Our estimates imply the Curie- Weiss constant to be 9 = 76.1 K. 

IV. CONCLUSIONS 

Several different series expansions have been calculated to high order for this model: 
high-temperature expansions, an Ising expansion at zero temperature, and dimer expansion 
at both zero and finite temperature. The first-order transition from the dimer phase to a 
Neel-ordered phase has been found to occur at J'/ J = 0.691(6), in good agreement with the 
original estimate 0.70(1) of Miyahara and Ueda @. The ground-state energy shows a sharp 
and distinct break in slope at that point, indicative of a first-order transition, and the triplet 
spin gap undergoes a small but definite discontinuity: the triplet spin gap is A = 0.14(5) at 
the transition point. A discontinuity in the Neel phase magnetization is less certain, but is 
not excluded by our results. 

The model has been fitted to the experimental susceptibility data for SrCu2 (1303)2 @, 
with parameters g = 2.108, J = 82.0 K and x = 0.678 if the energy gap is 19 K or g = 2.108, 
J = 83.2 K and x = 0.664 if the energy gap is 25 K. A detailed and accurate fit can be 
obtained at high temperature, and a reasonable fit at low temperatures is also obtained. 
The fitted value of x = (J'/J) is about 0.67, which is indeed very close to the transition 
point 0.691. 

We should summarize what is new in this work. The Ising expansion has been developed 
for this model for the first time, and this allows a more accurate determination of the phase 
transition point. We have also estimated the spontaneous magnetization in the Neel phase. 
The finite temperature dimer expansion has not been computed previously. This provides 
the most reliable theoretical estimate of the high temperature susceptibility and is used to 
estimate the experimental parameters of the real system. Apart from this we have also 
substantially increased the length of the dimer series at T = and the conventional high 
temperature expansion. The former has allowed us to present the first calculation of the full 
triplet excitation spectrum for this system. 
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FIGURES 



FIG. 1. (a) Lattice structure of the Cu 2+ spins of SrCu2(B03)2- The nearest-neighbor bonds 
are expressed by solid lines and the next-nearest-neighbor bonds by dashed lines, (b) Elemen- 
tary unit for interaction between a pair of nearest-neighbor bonds. The dotted lines denote the 
additional coupling ( J2) which still allow the perfect dimer state as eigenstate. 

FIG. 2. The ground-state energy per site Eq/NJ as function of J'/ J. The solid horizontal line 
is the energy of the perfect dimer state, while the solid points with error bars are the estimates 
from the Ising expansion. 

FIG. 3. The staggered magnetization M versus J/ J'. The solid points with error bars are the 
estimates from the Ising expansions. 

FIG. 4. Plot of triplet excitation spectrum A(k x , k y ) (derived from the dimer expansions) along 
high-symmetry cuts through the Brillouin zone for coupling ratios x = 0.6, 0.65, 0.678. 

FIG. 5. The triplet excitation gap A = A(0, 0) as a function of coupling ratios x derived from 
the dimer expansion. Several different integrated differential approximants to the series are shown. 
The errorbar indicates the gap at the critical point. 

FIG. 6. Internal energy per site E/NJ vs. temperature ksT '/ J for x = 1/2. The solid and 
dashed lines are the direct sum of the dimer expansion series at finite temperature to orders 
from 2 to 6, while the dotted lines are several different integrated differential approximants to the 
high-temperature series. 

FIG. 7. Specific heat per site C vs. temperature ksT/J for x = 0.4. The solid and dashed lines 
are the direct sum of the dimer expansion series at finite temperature to orders from 2 to 6, while 
the dotted lines are several different integrated differential approximants to the high-temperature 
series.. 

FIG. 8. Susceptibility \ per site x vs - temperature fceT/J for x = 1/2. The solid and dashed 
lines are direct sum of the dimer expansion series at finite temperature to orders from 2 to 6, while 
the dotted lines are several different integrated differential approximants to the high-temperature 
series.. 

FIG. 9. Internal energy per site E/NJ vs. temperature ksT/J for x = 3/2. The solid and 
dashed lines are direct sum of series of dimer expansion at finite temperature to order from 2 to 6, 
while the dotted lines are several different integrated differential approximants to high-temperature 
series. 
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FIG. 10. Susceptibility x P er site X vs - temperature ksT/J for x = 3/2. The solid and 
dashed lines are the direct sum of the dimer expansion series at finite temperature to orders 
from 2 to 6, while the dotted lines are several different integrated differential approximants to the 
high-temperature series.. 



FIG. 11. The two nearly horizontal solid lines in the lower window indicate where the energy 
gap is 19 K and 25 K, respectively. The other solid lines in both upper and lower windows indicate 
for a given value of J, the values of (x, g) which give a minimum value of P, and the corresponding 
value of P is also given in the upper window. The region bounded by two dashed lines in the 
lower window indicates the range in parameter space (J,x,g = 2.108) that P is less than or equal 
to 5P m i n , where P m i n is the minimum of P in the whole parameter space, indicated by the solid 
box point in the lower window. The dotted line in the lower window indicates where the Weiss 
constant 9 = (J + 4J')/4 = 92.5 K. 



FIG. 12. Comparison of the calculated temperature dependence of the susceptibility with ex- 
perimental data H (open points). The solid and dashed lines are direct sums of the dimer expansion 
series at finite temperature to orders from 2 to 6 for Heisenberg model with parameters J = 82 K, 
x = 0.678 and g = 2.108. 



FIG. 13. Same as Fig. 12, but the experimental data are compared with the estimates from 
high-temperature expansions (the solid lines representing various integrated differential approxi- 
mants to the high-temperature series.) 
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TABLES 



TABLE I. Series coefficients for the dimer expansion of the triplet excitation spectrum 
A(k x ,k y ) = JJ2k,n,m a k,n,m xk cos(mk x + nk y ). Nonzero coefficients ak, n ,m up to order k = 15 
are listed. 
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1 \ 


—2.169670151 x 10" 


-l 


/1 

(^10, 


i 


z 1 


6.075705255 x 10" 


-3 


( 4 


n 
U, 


n\ 
U J 


— 1 250000000x10" 


l 


/in 
(10, 


i 


i \ 

i y 


— 2 159589603x10" 


-4 


/ 1 /I 

(44, 


i 

1 




—9 902031964x10" 


-2 


/1 o 


2 


z 1 


5 358886719x10" 


-4 


( 5 


0, 


o) 


1 562500000x10" 


l 


/i 1 

(H, 


i 


i \ 

i y 


—6.345273723 x 10" 


-2 


/i r; 


i 

1 


1 \ 
-l y 


3.702751072 x 10" 


-2 


/I Q 

(^lo, 


n 
Z 


o^ 
z y 


1.970669593X 10" 


-3 


( 6 


U, 


U J 


2 343750000 X 10" 


2 


(12, 


i 


i \ 

i y 


— 1 767150368x10" 


-1 


/in 
(.4U, 


z 




— 7 866753472x10" 


-5 


/1 A 


n 
Z 


Z J 


3 433763327X 10" 


-3 


( 7 


0, 


0) 


— 3 687065072x10" 


1 


(13, 


1 


1) 


—2.169670151 x 10" 


-1 


(11, 


2 


-1) 


— 1 746396665x10" 


-4 


(15, 


2 


2) 


2.148257717x10" 


-3 


(8, 


o, 


0) 


-6.348922164x10" 


1 


(14, 


1 


1) 


-9.902031964x10" 


-2 


(12, 


2, 


-1) 


-3.807362686x10" 


-4 


(10, 


1 


-2) 


-7.866753472x10" 


-5 


(9, 


o, 


0) 


-3.904155213x10" 


1 


(15, 


1 


1) 


3.702751072x10" 


-2 


(13, 


2, 


-1) 


-1.039117707x10" 


-3 


(11, 


1 


-2) 


-1.746396665x10" 


-4 


(10, 


o, 


0) 


1.468893169x10" 


1 


(10, 


2, 


1) 


7.866753472x10" 


-5 


(14, 


2, 


-1) 


-2.816852119x10" 


-3 


(12, 


1 


-2) 


-3.807362686x10" 


-4 


(11, 


o, 


0) 


2.410994431x10" 


1 


(11, 


2, 


1) 


1.746396665x10" 


-4 


(15, 


2 


-1) 


-6.075705255x10" 


-3 


(13, 


1 


-2) 


-1.039117707x10" 


-3 


(12, 


o, 


0) 


-6.102249664x10" 


1 


(12, 


2 


1) 


3.807362686x10" 


-4 


(14, 


3, 


-1) 


-4.682591353x10" 


-7 


(14, 


1 


-2) 


-2.816852119x10" 


-3 


(13, 


o, 


0) 


-1.761368461 




(13, 


2, 


1) 


1.039117707x10" 


3 


(15, 


3, 


-1) 


-5.703426062x10" 


-7 


(15, 


1 


-2) 


-6.075705255x10" 


-3 


(14, 


o, 


0) 


-1.770029604 




(14, 


2, 


1) 


2.816852119x10" 


3 


(10, 


0, 


2) 


-8.680555556x10" 


-4 


(12, 


2 


-2) 


5.358886719x10" 


-4 


(15, 


o, 


0) 


-1.691084766x10" 


1 


(15, 


2, 


1) 


6.075705255x10" 


-3 


(11, 


0, 


2) 


-2.002495660x10" 


-3 


(13, 


2 


-2) 


1.970669593x10" 


-3 


(10, 


2, 


0) 


-2.486617477x10" 


4 


(14, 


3 


1) 


-4.682591353x10" 


-7 


(12, 





2) 


-7.904745973x10" 


-4 


(14, 


2 


-2) 


3.433763327x10" 


-3 


(11, 


2, 


0) 


-3.280526620x10" 


4 


(15, 


3 


1) 


-5.703426062x10" 


-7 


(13, 





2) 


6.061904493x10" 


-3 


(15, 


2 


-2) 


2.148257717x10" 


-3 


(12, 


2, 


0) 


9.612269774x10" 


4 


(6, 


1 


-1) 


-2.083333333x10" 


-2 


(14, 





2) 


1.576995064x10" 


-2 


(14, 


1 


3) 


-9.042245370x10" 


-6 


(13, 


2, 


0) 


4.323950651x10" 


3 


(7, 


1 


-1) 


-3.819444444x10" 


-2 


(15, 


o, 


2) 


1.850923153x10" 


-2 


(15, 


1 


3) 


-2.671619241x10" 


-5 


(14, 


2, 


0) 


6.550286294x10" 


3 


(8, 


1 


-1) 


-3.602430556x10" 


-2 


(10, 


1 


2) 


7.866753472x10" 


-5 


(14, 


1 


-3) 


-9.042245370x10" 


-6 


(15, 


2, 


0) 


8.172418539x10" 


4 


(9, 


1 


-1) 


-7.375458140x10" 


-3 


(11, 


1 


2) 


1.746396665x10" 


-4 


(15, 


1 


-3) 


-2.671619241x10" 


-5 


(6, 


1, 


1) 


-2.083333333x10" 


2 


(10, 


1 


-1) 


-2.159589603x10" 


-4 


(12, 


1 


2) 


3.807362686x10" 


-4 













TABLE II. Series coefficients of the Ising expansion for the ground-state energy per site 
Eq/(NJ) and the staggered magnetization M, for J' j J = 0.7,0.8,1,2 and t = 0. Nonzero co- 
efficients of A* up to order i = 9 are listed. 



t 


J'/J = 0.7 




J'/J = 0.8 




J'/J = 1 




J'/J = 2 












l 


E /(JN) 











-2.2500000000x10" 


i 


-2.7500000000x10" 


-3.7500000000x10" 


l 


-8.7500000000x10" 


l 


1 


0.0000000000 




0.0000000000 




0.0000000000 




0.0000000000 




2 


-2.2272727273x10" 


i 


-2.2857142857x10" 


l 


-2.5000000000x10" 


i 


-4.0000000000x10" 


i 


3 


1.0123966942x10" 


i 


8.1632653061x10" 


2 


6.2500000000x10" 


2 


4.0000000000x10" 


2 


4 


-1.3065004904x10" 


i 


-7.4186951609x10" 


2 


-3.5416666667x10" 


2 


-6.1067821068x10" 


3 


5 


2.6804197345x10" 


i 


1.2604414182x10" 


1 


5.0052083333x10" 


2 


1.0890524707x10" 


2 


6 


-5.3530295140x10" 


i 


-1.9791154855x10" 


1 


-6.1429832176x10" 


2 


-1.2407230180x10" 


2 


7 


1.1672953875 




3.2781449128x10" 


1 


7.4708752894x10" 


2 


8.8405316048x10" 


3 


8 


-2.9587106052 




-6.5298672657x10" 


1 


-1.1569637564x10" 


1 


-1.0969498862x10" 


2 


9 


7.7545598044 




1.3549175235 




1.8737141791x10" 


1 


1.0540011561x10" 


2 








Magnetization M 











5.0000000000x10" 


i 


5.0000000000x10" 


1 


5.0000000000x10" 


1 


5.0000000000x10" 


1 


1 


0.0000000000 




0.0000000000 




0.0000000000 




0.0000000000 




2 


-4.0495867769x10" 


i 


-3.2653061224x10" 


1 


-2.5000000000x10" 


1 


-1.6000000000x10" 


1 


3 


3.6814425244x10" 


i 


2.3323615160x10" 


1 


1.2500000000x10" 


1 


3.2000000000x10" 


2 


4 


-1.3553013187 




-6.2496725406x10" 


1 


-2.4361111111x10" 


1 


-5.1163207419x10" 


2 


5 


3.4305010340 




1.2470397038 




3.5369444444x10" 


1 


3.4797358726x10" 


2 


6 


-9.1962657747 




-2.5644449411 




-5.5283436053x10" 


1 


-4.6108122016x10" 


2 


7 


2.5590774140 xlO 1 




5.4148695622 




8.5367555708x10" 


1 


4.0349371222x10" 


2 


8 


-7.8821305450X10 1 




-1.3000427014X10 1 




-1.5762664485 




-5.9375271517x10" 


2 


9 


2.4067476616 xlO 2 




3.1003113436X10 1 




2.8719797054 




6.2701777926x10" 


2 
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TABLE III. Series coefficients for the dimer expansion at finite temperature of the magnetic 

(n) 

susceptibility x an d the logarithm of the partition function In Z. Nonzero coefficients ci / in Eq. 
(HH) up to order n = 6 are listed. 





,k,l) 


A; Z 


(n,k,l) 


k I 




(n 


k,l) 


k I 


(n,k,l) 


k I 




















x,f3) 














( °> 


0. 


0) 


4 


( 4 . 


3, 


1) 


-20689920 




( 5, 


2, 


3) 


-6519859200 


( 6, 


5, 


2) 


126045573120 


( o, 


0. 


1) 


-4 


( 4, 


4, 


1) 


-3239424 




( 5 . 


3, 


3) 


28138690560 


( 6, 


6. 


2) 


5370983424 


( 1, 


1. 


0) 


-64 


( 4 . 


1, 


2) 


23417856 




( 5 


4, 


3) 


1845657600 


( 6, 


1, 


3) 


-386640691200 


( !> 


1. 


1) 


128 


( 4 . 


2 : 


2) 


50377728 




( 5 


5, 


3) 


—56417280 


( 6, 


2, 


3) 


-2967429427200 


( 1, 


1. 


2) 


-64 


( 4, 


3, 


2) 


56881152 




( 5 . 


1, 


4) 


-1577410560 


( 6, 


3, 


3) 


-1905248194560 


( 2 > 


1. 


0) 


640 


( 4, 


4, 


2) 


1916160 




( 5 . 


2, 


4) 


6189373440 


( 6, 


4, 


3) 


-2606694359040 


( 2 > 


2, 


0) 


1344 


( 4 . 


1, 


3) 


-22579200 




( 5 


3, 


4) 


-36834969600 


( 6, 


5, 


3) 


9843886080 


( 2 ; 


1. 


1) 


-4224 


( 4 . 


2, 


3) 


-109221888 




( 5 


4, 


4 ) 


1241210880 


( 6, 


6, 


3) 


-5610921984 


( 2 > 


2, 


1) 


-3456 


( 4, 


3, 


3) 


-68871168 




( 5. 


5, 


4) 


33269760 


( 6, 


1, 


4) 


412034826240 


( 2 ! 


1, 


2) 


7680 


( 4, 


4, 


3) 


1858560 




( 5. 


1, 


5) 


188743680 


( 6, 


2, 


4) 


3178633052160 


( 2 ! 


2, 


2) 


2880 


( 4 . 


1, 


4) 


7593984 




( 5 


2, 


5) 


-2117468160 


( 6, 


3, 


4) 


4207696773120 


( 2 > 


1, 


3) 


-4096 


( 4, 


2. 


4) 


105916416 




( 5 . 


3, 


5) 


23898193920 


( 6, 


4, 


4) 


3593735331840 


( 2 ! 


2, 


3) 


-768 


( 4, 


3, 


4) 


37969920 




( 5 . 


4, 


5) 


-2151383040 


( 6, 


5, 


4) 


-331372339200 


( 3 > 


1. 


0) 


16128 


( 4 . 


4, 


4) 


-2433024 




( 5 


5, 


5) 


36126720 


( 6, 


6, 


4) 


15752116224 


( 3 ! 


2, 


0) 


-39168 


( 4 . 


2 : 


5) 


—37355520 




( 5 


3, 


6) 


-6102712320 


( 6, 


1, 


5) 


— 160650362880 


( 3, 


3, 


0) 


-36864 


( 4, 


3, 


5) 


-7667712 




( 5, 


4, 


6) 


747110400 


( 6, 


2, 


5) 


-1772148326400 


( 3, 


1. 


1) 


-103680 


( 4, 


4, 


5) 


712704 




( 5, 


5, 


6) 


-21626880 


( 6, 


3, 


5) 


-5437322035200 


( 3 > 


2, 


1) 


359424 


( 5, 


1, 


°) 


-190402560 




( 6 


1, 


°) 


-3039897600 


( 6, 


4, 


5) 


-2693322178560 


( 3 > 


3, 


1) 


99072 


( 5, 


2, 


0) 


-75018240 




( 6, 


2, 


0) 


35239587840 


( 6, 


5, 


5) 


441996410880 


( 3 ) 


1, 


2) 


179712 


( 5. 


3, 


o) 


-157562880 




( 6, 


3, 


0) 


24536862720 


( 6, 


6, 


5) 


-22042312704 


( 3, 


2, 


2) 


-933120 


( 5. 


4, 




-142878720 




( 6 


4, 


0) 


17376768000 


( 6, 


1, 


6) 


9059696640 


( 3, 


3, 


2) 


-78336 


( 5, 


5, 


0) 


-43966464 




( 6 


5, 


0) 


8583413760 


( 6, 


2, 


6) 


402165596160 


( 3, 


1, 


3) 


-92160 


( 5, 


1. 


1) 


1436820480 




( 6, 


6, 


0) 


1929157632 


( 6, 


3, 


6) 


3619867852800 


( 3, 


2, 


3) 


944640 


( 5, 


2, 


1) 


-283760640 




( 6 


1, 


1) 


-2405099520 


( 6, 


4, 


6) 


1007563898880 


( 3, 


3, 


3) 


6912 


( 5, 


3, 


1) 


2067609600 




( 6 


2, 


1) 


-386275184640 


( 6, 


5, 


6) 


-240553820160 


( 3, 


2. 


4) 


-331776 


( 5, 


4, 


1) 


1148129280 




( 6 


3, 


1) 


-162191831040 


( 6, 


6. 


6) 


13220904960 


( 3, 


3, 


4) 


9216 


( 5, 


5, 


1) 


119927808 




( 6 


4, 


1) 


-204870942720 


( 6, 


3, 


7) 


-958440407040 


( 4 , 


1. 


0) 


1290240 


( 5, 


1, 


2) 


-3521018880 




( 6 


5, 


1) 


-63168215040 


( 6, 


4, 


7) 


-140236554240 


( 4 , 


2. 


0) 


1256448 


( 5, 


2, 


2) 


2806732800 




( 6 


6, 


1) 


-5717993472 


( 6, 


5, 


7) 


48625090560 


( 4 , 


3. 


0) 


2377728 


( 5, 


3, 


2) 


-11009249280 




( 6, 


1, 


2) 


131641528320 


( 6, 


6, 


7) 


-2901934080 


( 4 , 


4, 


0) 


1185024 


( 5, 


4, 


2) 


-2687846400 




( 6, 


2, 


2) 


1509814702080 










( 4 , 


1. 


1) 


-9722880 


( 5, 


5, 


2) 


-67313664 




( 6 


3, 


2) 


611100979200 










C 4 ) 


2. 


1) 


-10973184 


( 5 . 


1, 


3) 


3663267840 




( 6 


4, 


2) 


1026448035840 


























InZ 


X, 


3) 














( 2. 


1. 


0) 


-576 


( 4 . 


2, 


2) 


-21150720 




{ 5 


5, 


2) 


-16692480 


( 6, 


4, 


2) 


151663104000 


( 2 > 


2, 


0) 


576 


( 4 . 


3, 


2) 


10824192 




( 5 


1, 


3) 


-313528320 


( 6, 


5, 


2) 


-29455488000 


( 2 ! 


1. 


1) 


2880 


( 4, 


4, 


2) 


-1534464 




( 5. 


2, 


3) 


1816888320 


( 6, 


6. 


2) 


2645913600 


( 2 i 


2, 


1) 


-1152 


( 4 . 


1, 


3) 


2322432 




( 5 


3, 


3) 


58475520 


( 6, 


1, 


3) 


143954288640 


( 2 > 


1, 


2) 


-2304 


( 4 . 


2, 


3) 


22975488 




( 5 


4, 


3) 


-145981440 


( 6, 


2, 


3) 


262998005760 


( 2 > 


2, 


2) 


576 


( 4, 


3, 


3) 


-12441600 




( 5 . 


5, 


3) 


15068160 


( 6, 


3, 


3) 


676278097920 


( 3 ! 


1, 


0) 


-10368 


( 4, 


4, 


3) 


1465344 




( 5 . 


2, 


4) 


-663552000 


( 6, 


4, 


3) 


-444080517120 


( 3 J 


2, 


0) 


-10368 


( 4 . 


2, 


4) 


-9289728 




( 5 


4, 


4) 


51425280 


( 6, 


5, 


3) 


100333209600 


( 3 1 


3, 


0) 


1728 


( 4 . 


3, 


4) 


4644864 




( 5 


5, 


4) 


-4700160 


( 6, 


6, 


3) 


-8033458176 


( 3 ! 


1. 


1) 


51840 


( 4, 


4, 


4) 


-470016 




( 6, 


1, 


0) 


-6057815040 


( 6, 


1. 


4) 


-53031075840 


( 3, 


2, 


1) 


20736 


( 5. 


1, 


0) 


87713280 




( 6 


2, 


0) 


-5008988160 


( 6, 


2, 


4) 


-219410104320 


( 3, 


3, 


1) 


-3456 


( 5, 


2 : 


°) 


-42716160 




( 6 


3, 


°) 


-8885790720 


( 6, 


3, 


4) 


-814218117120 


( 3, 


1, 


2) 


-41472 


( 5, 


3, 


0) 


-54743040 




( 6, 


4, 


0) 


-1159557120 


( 6, 


4, 


4) 


620258549760 


( 3, 


2, 


2) 


-10368 


( 5, 


4, 


0) 


-2073600 




( 6, 


5, 


0) 


187038720 


( 6, 


5, 


4 ) 


-143287418880 


( 3, 


3, 


2) 


1728 


( 5, 


5, 


0) 


-656640 




( 6 


6, 


0) 


-37324800 


( 6, 


6, 


4) 


10433193984 


( 4, 


1. 


°) 


-705024 


( 5, 


1. 


1) 


-516948480 




( 6 


1, 


1) 


49705436160 


( 6, 


2. 


5) 


66886041600 


( 4, 


2, 


0) 


-1575936 


( 5, 


2. 


1) 


664796160 




( 6, 


2, 


1) 


38078760960 


( 6, 


3, 


5) 


537954877440 


( 4, 


3, 


0) 


-124416 


( 5, 


3, 


1) 


167961600 




( 6, 


3, 


1) 


90482780160 


( 6, 


4, 


5) 


-419908976640 


( 4, 


4, 


0) 


-13824 


( 5, 


4, 


1) 


-38983680 




( 6 


4, 


1) 


-17612328960 


( 6, 


5, 


5) 


94038589440 


( 4, 


1. 


1) 


4105728 


( 5, 


5, 


1) 


6981120 




( 6, 


5, 


1) 


1647267840 


( 6, 


6, 


5) 


-6303744000 


( 4, 


2. 


1) 


9040896 


( 5, 


1, 


2) 


742763520 




( 6, 


6, 


1) 


-164395008 


( 6, 


3, 


6) 


-147786301440 


( 4, 


3, 


1) 


-2903040 


( 5, 


2, 


2) 


-1775416320 




( 6 


1, 


2) 


-134570833920 


( 6, 


4, 


6) 


110839726080 


( 4, 


4, 


1) 


552960 


( 5, 


3, 


2) 


-171694080 




( 6 


2, 


2) 


-143543715840 


( 6, 


5, 


6) 


-23463198720 


( 4, 


I, 


2) 


-5723136 


( 5, 


4, 


2) 


135613440 




( 6 . 


3, 


2) 


-333825546240 


( 6, 


6. 


6) 


1459814400 



14 



